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Abstract. A 7-ray line production calculation in astro- 
physics depends on i) the composition and energy source 
spectrum of the energetic particles, ii) the propagation 
model, and iii) the nuclear cross sections. The main diffi- 
culty for model predictions and data interpretation comes 
from the fact that the spectrum of the particles which ac- 
tually interact in the ISM - the propagated spectrum, is 
not the same as the source spectrum coming out of the ac- 
celeration site, due to energy-dependent energy losses and 
nuclear destruction. We present a different approach to 
calculate 7-ray line emission, based on the computation 
of the total number of photons produced by individual 
energetic nuclei injected in the interstellar medium at a 
given energy. These photon yields take into account all 
the propagation effects once and for all, and allow one to 
calculate quickly the 7-ray line emission induced by en- 
ergetic particles in any astrophysical situation by using 
directly their source spectrum. Indeed, the same photon 
yields can be used for any source spectrum and composi- 
tion, as well as any target composition. In addition, these 
photon yields provide visual, intuitive tools for 7-ray line 
phenomenology. 

Key words: Cosmic rays; Gamma rays: theory; Nuclear 
reactions 



1. Introduction 

High energy astrophysics is experiencing a considerable 
development, notably through the operation of a new gen- 
eration of gamma-ray observatories, both on the ground 
and onboard satellites. The analysis of gamma-ray emis- 
sion from compact and diffuse sources is one of the most 
valuable ways to study high energy processes in the uni- 
verse, and it is expected to put stronger and stronger 
constraints on the models in the near future. Among the 
processes of interest, the production of gamma-ray lines 
by energetic particles (EPs) interacting in the interstel- 
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lar medium (ISM) has received increased attention in the 
last few years, in connection with the data of the Comp- 
ton Gamma-Ray Observatory as well as the forthcoming 
satellite INTEGRAL. In addition to the study of the high 
energy sources, the study of EP interactions is important 
for the understanding of the EPs themselves, of which the 
Galactic cosmic rays (GCRs) are one of the main compo- 
nents. 

Information about the acceleration sites and processes 
is also provided by the determination of the EP energy 
spectrum and chemical composition, which can be derived, 
in principle, from the measurement of gamma-ray line ra- 
tios and profiles. However, the information contained in 
the gamma-ray observational data relates to the spectrum 
and composition of the propagated particles (i.e. those who 
actually interact in the ISM), not the source particles. The 
difference arises from the fact that the EPs experience var- 
ious types of interactions while they 'propagate' from their 
source to the place where they produce gamma-ray lines. 
In particular, they lose energy through Coulombian in- 
teractions in a way which depends on both their energy 
and chemical nature, so that the propagated population 
of EPs is not identical to the source population, freshly 
coming out of the acceleration process. 

Ideally, one would divide the whole process into three 
successive stages (e.g. Parizot and Lehoucq, 1999): par- 
ticle acceleration, propagation and interaction, with the 
gamma-ray production arising during the last stage only. 
In reality, of course, the reactions leading to gamma-ray 
production occur all the time, from the injection of the 
EPs into the ISM until they have slowed down to energies 
below the interaction thresholds. It is therefore necessary 
to sum the contributions of all the instants following ac- 
celeration (the gamma-ray emission occurring during ac- 
celeration itself can usually be neglected, except for the 
rarest, highest energy particles, which spent a long time in 
the accelerator). In steady state situations, this is equiva- 
lent to calculating the equilibrium distribution of EPs and 
integrating the relevant cross sections over this so-called 
propagated distribution. 
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From a technical point of view, the most difficult part 
consists in calculating the propagated spectrum, taking 
into account the energy losses and the energy-dependent 
escape time of the particles out of the confinement re- 
gion, where the gamma-ray production is evaluated. This 
is what prevents a straightforward calculation of the ex- 
pected gamma-ray line fluxes from the knowledge of the 
nuclear cross sections and the EP source distribution. 
Therefore, we propose here to work out this step once 
and for all, in the case of a steady state and a thick tar- 
get, by calculating the integrated effect of energy losses 
on individual particles injected at any energy in the ISM. 
In Sect. ||, we shall justify the fact that the 'propagation 
step' in a standard 7-ray line emission calculation can be 
'factorized out' and calculated separately, independently 
of the EP source spectrum and composition. This is true 
if the mctallicity of the propagation medium does not ex- 
ceed several tens of times the solar metallicity, as in most 
of the astrophysically sensible situations. As a result, we 
shall obtain the absolute 7-ray yields of energetic nuclei 
as a function of their initial energy, from which the 7-ray 
line emission induced by EPs of any spectrum and compo- 
sition can be straightforwardly calculated. In addition to 
making 7-ray line calculations much easier, these absolute 
yields (or particle efficiencies for 7-ray line production) 
considerably help phcnomcnological interpretation of the 
observational data, as these yields only need to be con- 
volved with the EP source spectra, rather than propagated 
ones. 

2. Standard gamma-ray line calculations 

In order to calculate the gamma-ray line emission in a 
given region of the ISM, one needs to integrate the nu- 
clear excitation cross sections over the local flux of en- 
ergetic particles, and sum all the contributions to each 
gamma-ray line. If i represents the projectile, j the tar- 
get nucleus, and k the excited nucleus produced in the 
interaction between i and j, or equivalently the 'photon 
species' emitted through nuclear de-excitation, the 7-ray 
line emission rate, in ph/cm 3 /s, from all the i + j — > k 
nuclear reactions, reads: 

dTV*. ^ f + °° 

= J Ni{E) [ nj a itj , k (E)v{E)} dE 

iJ \ (1) 
= V / $i{E)njO- i>:j;k (E)dE, 
id Jo 

where Ni(E) is the spectral density of the projectiles i, 
in cm" 3 (MeV/n)" 1 , <&i(E) = Ni(E) x v(E) is the corre- 
sponding flux, in cm" ^s-^MeV/n)" 1 , v(E) is the velocity 
of the projectiles (independent of the nuclear species, i, if 
the energy E is expressed in MeV/n), rij is the number 
density of the target nuclear species j, and Cij-k is the 
cross section for the reaction i + j — > k. For instance, k 



might represent photons from the 12 C de-excitation line 
at 4.44 MeV, produced by the reaction p + 12 C — ► 12 C* 
ora + 16 — > 12 C*. 

Assuming that the nuclear excitation cross-sections are 
known, the main challenge is to estimate the EP fluxes, 
for each nuclear species, in the gamma-ray source. This 
depends on the acceleration process at the origin of the 
injection of EPs in the ISM, and on what happens to the 
particles once they leave the accelerator, which involves 
the energy losses, the rate of escape from the region con- 
sidered, and the particle destruction in inelastic nuclear 
processes. As far as the acceleration is concerned, it can 
be characterized here by the so-called injection function, 
Qi(E), which gives the number of particles of species i 
injected at energy E in the ISM (i.e. leaving the accel- 
eration process and not being further accelerated after- 
wards), in cm _3 s -1 (MeV/:n) _1 . The function Qi(E) will 
either be taken as the outcome of some particular accel- 
eration model (e.g. diffusive shock acceleration), or phe- 
nomenologically postulated so as to reproduce some par- 
ticular observation (e.g. from INTEGRAL data). 

Most studies so far have assumed that the shape of the 
injection spectrum, Qi(E), was independent of the nuclear 
species, and could be re- written as aiQ(E), where the iso- 
topic abundances a.i characterize the chemical composi- 
tion of the EPs at the source. However, this simplification 
is not required and one will allow here for a different spec- 
trum for each nuclear species, which is equivalent to an 
energy dependent EP composition. 

In order to use Eq. (|]J) to calculate the gamma-ray 
emission produced in the region under consideration, one 
needs to derive the EP fluxes, $i(E) or N{(E), from the 
injection functions, Qi(E), supposed known. The standard 
way to do this has been described in Parizot and Lehoucq 
(1999, and references therein) for the general case where 
the injection function as well as the conditions of propaga- 
tion are time-dependent. It consists in solving the so-called 
propagation equation, which takes the following form in 
the stationary case: 

± {Ei (E)N i (E)) = Q i (E)-^- y (2) 

Here, Ei(E) is the energy loss function, in 
(MeV/n)s _1 , giving the rate of energy loss for nuclei 
of species i in the considered propagation medium, and 
T^ ot (E) is the total 'loss time' taking into account nuclear 
destruction and particle escape out of the region consid- 
ered. It can be expressed in terms of the total inelastic 
cross sections for nuclei of species i in a medium made 
of j nuclei alone, af^ st (E), and the mean escape time, 
rf c (£), as: 

^ = ^+Y,^t est ( E >w> ( 3 ) 
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Fig. 1. Energy losses per unit grammage of matter passed through, in (MeV/n)/(g cm -2 ): on the left, for a 12 C nucleus 
in various chemically pure propagation media indicated by the labels; on the right, for various nuclei indicated by the 
labels in a propagation medium made of pure Hydrogen. 



where rij is as above the number density of nuclei of species 
j in the propagation medium. 

In Eq. (||), the injection function Qi{E) acts as a source 
term, while the equilibrium EP distribution iVj (E) is what 
we want to calculate. The formal solution of the stationary 
propagation equation reads: 



\Ei(E)\ Je 



Qi(E ia )V i (E in ,E)dE il 



(4) 



where Vi{E- ln , E) can be interpreted as the survival proba- 
bility, in the propagation medium considered, of a particle 
injected at energy E- m and losing energy down to energy 
E. It obviously depends on the total loss time at each 
energy between E m and E and the energy loss function, 
and can be expressed as follows (see Parizot and Lehoucq, 
1999): 



(5) 



3. Approached EP propagation universality 

Equation ([!]) allows one to calculate the photon emis- 
sion rate for any gamma-ray de-excitation line for which 
the production cross sections are known. It is based on 
the computation of the EP distribution function given by 
Eq. (^) , which in turn requires the knowledge of the energy 
losses and the escape and destruction times. 



3.1. Energy loss rates 

Let us first consider the energy loss rate, Ei(E), for par- 
ticles of species i. It can be expressed as a sum over the 
nuclear species present in the propagation medium: 



E i {E) = Y,n j 



AE 
dx 



v(E), 



(6) 



where (dE/dx)j is the energy loss per unit grammage in 
a medium of pure j nuclei, in (MeV/n)/(gcm -2 ). In this 
expression, (dE/dx)j is a 'universal function' of energy 
which can ideally be calculated from first physical prin- 
ciples or, failing that, extrapolated from laboratory mea- 
surements, while the astrophysics comes in the chemical 
composition, the degree of ionization and the density of 
the propagation medium, rij . Some examples of energy loss 
functions used in this paper are shown in Fig. |l[ They have 
been calculated using the program of J. Kiener (1994), 
based on the modified Bethe formula taking into account 
the effective charge of the projectile. This program imple- 
ments the Ziegler tables for the various stopping powers, 
as corrected by Hubert et al. (1989) according to a semi- 
empirical procedure, where a new parameterization for the 
effective charge is deduced from a very large set of exper- 
imental stopping power values in the range 3-80 MeV/n. 
The typical error in the energy loss function is between 2% 
and 10%, i.e. less than the uncertainty about the nuclear 
excitation cross-sections. 

The energy range of interest for the calculation of 
gamma-ray line emission is between a few McV/n (cor- 
responding to the nuclear excitation thresholds) and a 
few hundreds of MeV/n. At higher energy, the contribu- 
tion of the EPs to the gamma-ray line emission is small, 
because of their reduced number (decreasing power-law 
source spectrum) as well as because they are destroyed by 
nuclear reactions before they reach the peaks of the nu- 
clear excitation cross sections, as further discussed below. 

In all the calculations presented here, the ambient 
medium is assumed neutral, which may not be appropri- 
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Fig. 2. Stationary energy spectra of 1 H (left) and 16 (right) nuclei injected with a simple power-law spectrum in 
energy of index 2, after propagation in different media. The spectra corresponding to propagation media with a 
metallicity up to 10 Zn are virtually indistinguishable from the propagated spectrum in a metal-free gas. 



ate for a number of astrophysical situations. However, one 
can estimate that the effect of an ambient ionized medium 
is small, except for the lowest energies. The reason why 
the energy losses depend on the ionization state of the 
propagation medium is that it is more difficult for an en- 
ergetic ion to capture a free electron than to capture an 
electron from an atom at rest. Indeed, in the latter case, 
the orbital motion of the electron reduces the velocity dif- 
ference between the energetic ion and the electron, and 
thereby facilitates capture. As a consequence, the equilib- 
rium between electron stripping and electron capture de- 
pends on the ambient medium, and an energetic atom is 
on average more ionized when it travels through a plasma 
than through a neutral medium (Chabot et al., 1995a). 
This results in a higher effective charge, and thus a higher 
stopping power (or larger energy losses). However, when 
the projectile is too energetic, its relative velocity with 
even orbital electrons is too high for charge exchange to 
be efficient anyway. Therefore, the difference between an 
ionized and a neutral ambient medium becomes negligible, 
and the energy losses are almost identical. From the quan- 
titative point of view, the stopping power of a plasma is 
higher than that of a neutral medium by a factor of about 
40 below 100 keV/n, but only 2 or 3 at a few MeV/n, and 
their difference is negligible above 100 MeV/n (Hoffman et 
al., 1994; Chabot et al., 1995b). Therefore, we will assume 
throughout that the propagation medium is neutral. 

As far as the chemical composition is concerned, 
in most astrophysically relevant cases the propagation 
medium will be but the ISM, whose composition is rel- 
atively well known (Anders and Grevesse, 1989). Now, al- 
though the heavy elements are of course crucial to the 
calculation of gamma-ray de-excitation line emission, the 
ISM is so much dominated by H and He nuclei that one can 
neglect all other elements in Eq. (||) , and calculate the en- 



ergy loss rate as if the ISM were made simply of 91% of H 
and 9% of He (by number). To demonstrate this, we have 
calculated the propagated (equilibrium) spectrum of the 
different energetic nuclei subject to energy losses in me- 
dia of various metallicity. Results are shown in Fig. ^ for 
metallicities ranging from (H and He only) to 1000 times 
solar. A significant change in the particle propagated spec- 
trum can only be noticed for ambient metallicities larger 
than ten times the solar metallicity, Zq. Since most of the 
astrophysically relevant media are not that rich in met- 
als, we will assume that the propagation of the EPs is 
independent of metallicity. Note that even pure SN ejecta 
have a metallicity less than 30 times Zq, so that even in 
such a metal-rich medium, neglecting the interaction of 
the EPs with the metals as they propagate through the 
ambient medium will lead to an error smaller than 20% 
on the propagated particle distribution, Ni{E), and thus 
also on the gamma-ray line production rates. 

3.2. Total inelastic cross sections and survival 
probabilities 

The second crucial physical ingredient necessary to calcu- 
late the EP propagated spectra is the total 'catastrophic 
loss time', T* ot (i?), including both nuclear destruction and 
escape (see Eq. (||)). The latter will be neglected here, 
which amounts to say that we consider a distribution of 
EPs interacting with a thick target. This is well justified 
for the relatively low energy particles which are respon- 
sible for most of the gamma-ray line emission, as their 
range is of the order of 1 g/cm 2 , i.e. significantly less than 
the typical escape grammage of cosmic rays. One can thus 
replace the total catastrophic loss timescale, r* ot (i?), by 
the destruction timescale, Tf(E). 
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Fig. 3. EP survival probability down to 10, 30, 50 and 
100 MeV/n, as a function of the injection energy, E m . 

Concerning the EP destruction through nuclear reac- 
tions, we use semi-empirical total inelastic cross sections 
from Silberberg et Tsao (1990). As for the energy loss 
function, we have just seen that the EP destruction time 
can be calculated as if the propagation medium were made 
of pure H and He, to an excellent approximation (even for 
relatively high ambient metallicities) . In all of our calcula- 
tions, the main uncertainty comes from the relatively poor 
knowledge of the nuclear cross sections, rather than from 
the use of a 'universal' (metal-free) propagation medium. 

The joint knowledge of E^E) and T^ ot (E) allows one 
to calculate the survival probability of an EP injected at 
any energy Ei n while it loses energy down to the energy E 
at which it interacts with the ISM to produce a gamma- 
ray. It can be useful to define the energy loss timescale, 
t} oss {E) for an EP of type i as: 



ni° ss (E) = E/\Ek(E)\ 



(7) 



The survival probability given by Eq. (|J) can then be 
rewritten as: 



Vt(E ta ,E) = exp[ - 



t} oss {E') dE' 



(8) 



Fig. H shows the results obtained in the case of a 12 C 
nucleus. As can be seen, particles injected at energies 
greater than 1 GeV/n have very little chance to survive 
long enough to reach energies of a few tens of MeV/n, 
where the nuclear excitation cross sections are maximum. 
As a consequence, their contribution to the total gamma- 
ray emission rate will be negligible, not mentioning the 
fact that higher energy particles are usually less numerous 
than lower energy ones, for the most common, power-law 
source spectra. This effect reminds us that nuclear exci- 
tation is only one of a number of competing nuclear pro- 
cesses affecting an energetic particle, so that the efficiency 



of gamma-ray line production depends on the balance be- 
tween several cross sections. 



4. Gamma-ray yields of individual EPs 

We now make use of the results of the previous section to 
derive an alternate way to calculate gamma-ray line pro- 
duction from EPs. This idea is the following: since particle 
propagation is 'universal' (i.e. independent of the physical 
conditions of the propagation medium, as long as metals 
can be neglected), there should be a way to work it out 
once and for all. This way is the following: instead of cal- 
culating the equilibrium EP distribution function (after 
propagation) , and then integrate the nuclear cross section 
over this distribution, one can calculate the total pho- 
ton yield, from injection to rest, of one EP of a given 
kind thrown in the ISM at any given initial energy, inde- 
pendently of the other particles (i.e. independently of the 
global spectrum and composition of the EPs), and then 
integrate the individual photon yields over the source dis- 
tribution function of each kind of EPs accelerated. The im- 
portant point is that the integration is now over the source 
distribution function, which is known from acceleration 
models or postulated in a phenomenological study, rather 
than over the propagated distribution function, which no 
longer has to be calculated. 

From the mathematical point of view, the above idea 
amounts to a simple change of the order of two integra- 
tions. Indeed, combining Eqs. (|l|) and (|J), one can rewrite 
the 7-ray emission rate as follows (specializing to one nu- 
clear reaction for illustration): 



dA 7 

~dT 



= / dE 



dE, 



n a(E)v(E) 
\E(E)\ 



Q(E in )T(E in ,E), (9) 



where hq is the density of the propagation medium and 
Q{E) is the EP source spectrum. Now, as is made obvious 
by Fig. this expression can be re-written as: 



d7V 7 
~dT 



+ dE n ( E dE no < E)v{E) Q{E- m )V{E in ,E). 



\E{E)\ 



(10) 



Getting the source function, Q(Ei n ), out of the integral 
over E, one obtains the following expression for the 7-ray 
emission rate (adding the contribution of all the reactions 
involved): 



dA 7 
~df 



r+00 

V / Q l {E ia )a 3 N 1 , rn {E- m )dE- ln , (11) 



where aj — rij/riQ is the relative abundance of nuclei of 
species j in the target, and 



■M j'n (-^in) 



Eia n a ijn (E)v(E) 
\E(E)\ 



Vi(E in ,E)dE. (12) 
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Fig. 5. On the left (a): Absolute photon yields, A/i J;7 , in the 6.13 MeV line of 16 through various channels, as a 
function of the injection energy of the projectile. The latter is the first nucleus appearing in the label, and the target 
is the second. On the right (b): Gamma-ray production efficiency, in photon/erg, for the 6.13 MeV line of 16 through 
various channels, as a function of the injection energy of the projectile. 
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Fig. 4. Graphical demonstration of the equivalence be- 
tween Eq. (^|) and Eq. (|io|): the shaded area is the inte- 
gration domain, divided into vertical and horizontal slices, 
respectively. 

The physical interpretation of Mi^iEya) is straight- 
forward: it is the number of photons of species 7 produced 



in a target made solely of nuclei of species j, by one pro- 
jectile of species i injected in the ISM at the energy E- mi 
integrated over its entire life (i.e. from its injection un- 
til it has lost so much energy that it is below the nuclear 
excitation threshold). Note that the lower bound of the in- 
tegral can be replace by the energy threshold of the cross 
sections. Now the interesting point is that the absolute 
photon yields, Ni^- n {E- m ), can be calculated from physi- 
cal quantities alone and is independent of astrophysics: as 
can be seen from Eqs. dl2| ) and (jjj), it only depends on the 
nuclear cross sections and energy loss rates. These can be 
calculated or measured once and for all, and so is it for 

As anticipated, the great advantage of this formulation 
is that once the quantities A/i.j ;7 (£ , i n ) have been calcu- 
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Fig. 7. Same as Fig. || for the 12 C line at 4.44 MeV, for reactions involving H nuclei (top) and He nuclei (bottom). 



lated, the actual 7-ray emission rate in a given astrophys- 
ical situation can be derived from Eq. (|ll]) which gathers 
all the astrophysical information (namely the EP spec- 
trum and composition, and the target composition), but 
which is now expressed in terms of the source spectrum, 
rather than the propagated one. To better understand the 
signification of this transformation, it suffices to compare 
Eqs. (|l|) and (^lj). We have replaced the propagated spec- 
tral density of the EPs, Ni(E), by their injection function, 
Qi(E), and the cross sections o~ij-j by our absolute pho- 
ton yields, A/i J;7 , which play the role of 'effective cross- 
sections' (although their physical dimension is different) 
taking into account the propagation of the EPs in the am- 
bient medium. It should be stressed that the individual 
photon yields behave as universal physical quantities and 
can be used with any source spectrum, any EP composi- 
tion and any target composition. 

Two comments are in order here. First, the above ex- 
pression giving the photon yields Ni^- tl {E- ln ) may seem 
to depend on the density, no, of the propagation medium 



(e.g. the ISM). This is actually not the case, as the energy 
loss rate appearing in the denominator is also proportional 
to this density. The second comment concerns the univer- 
sality of EP propagation, which is crucial in the approach 
developed here. Indeed, in principle the EP energy loss 
rates and survival probabilities depend on the propaga- 
tion medium, so that a different photon yield should be 
calculated for each propagation medium. These specific 
photon yields could still be used with any source spectrum 
and composition, but not with any target composition, as 
the latter is usually the same as that of the propagation 
medium. However, as shown in Sect. ^, the dependence of 
the energy loss rates and survival probabilities with metal- 
licity is negligible in most situations, so that the photon 
yields A/i J;7 can indeed be considered universal. 



5. Results and emission rates reconstruction 
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Fig. 8. Same as Figs. §for the 20 Ne line at 1.63 MeV and the 56 Fe line at 0.847 MeV. 



5.1. Nuclear excitation cross sections 

Before showing the individual photon yields calculated for 
the main expected gamma-ray lines in the ISM, we should 
say a word about cross sections, as they provide the major 
uncertainty in our results and most of them are extrap- 
olated from relatively scarce experimental data. We have 
used the data from Ramaty et al. (1979), updated with 
more recent experimental data whenever possible (Dyer et 
al., 1981,1985; Lang et al., 1987; Lesko et al. 1988; Kiener 
et al., 2001). The paper by Kiener et al. (2001, and pri- 
vate communication) has been used to derive a general 
procedure to extrapolate the data at high energy. 

Above the resonance peak of the cross sections, when- 
ever the data is missing we assume that the cross section 
obeys a simple law in a x E~ x + b (above 20-30 MeV/n, 
say). The best fit of the data for the 12 C(p,p7) reac- 
tion gives o(E) = 24337 x E~ 1J4 + 4 mbarn, where E 
is in MeV. For the 12 C(a, aj) reaction, one finds cr(E) = 
31400 x _E -1 - 51 + 7 mb, and for the excitative spallation re- 



action 16 0(p, pa) 12 C|. 44 , a{E) = 75951 x S" 1 - 98 + 4.3 mb 
(Kiener, 2001; private communication). 

As a first approximation, we assume that the excita- 
tion cross sections of other nuclei (namely, 14 N, 16 0, 20 Ne 
and 56 Fe) obey the similar laws with the same value of 
the exponent x, i.e. 1.74 and 1.51 respectively for proton- 
induced and alpha-particle-induced excitations. The con- 
stant value at high energy, B, is simply taken as the value 
of B for the 12 C excitation cross sections, but scaled to 
the measured value of the cross section at the resonance 
peak (i.e. proportionally to the peak value). Finally, the 
value of A is obtained by imposing the continuity of the 
cross sections above the peak. The same procedure is ap- 
plied to excitative spallation reactions (i.e. an exponent of 
1.98 and a high energy value scaled proportionally to the 
peak value), although such an extrapolation is more prob- 
lematic in this the 16 O nucleus has an atypical 
structure with a large component of four a particles. 

The above procedure can be summarized by the fol- 
lowing expression giving the cross section above the peak 
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Table 1. Data relative to the fit of the cross section as 
given by Eq. [l^. The reaction is given in the first column, 
with the names of the target and projectile and the energy 
of the resulting gamma ray in MeV. The second and third 
columns give the energy of the cross section peak and the 
corresponding value of the cross section. 



reaction 


-t-^pcak 


<^~pcak 


ITT i 5(jTTV, n 0/I7 

ri+ re={j.o4:( 


ion 


OD1 


4 Hp+ 56 Pp— Rd7 


-LO.U 




1 H+ 20 Ne=1.634 


11.0 


330 


4 He+ 20 Ne=1.634 


3.19 


232 


1 H+ 24 Mg=1.634 


23.0 


180 


1 H+ 28 Si=1.634 


21.7 


157 


1 H+ 14 N=2.313 


10.0 


106 


4 He+ 14 N=2.313 


3.25 


119 


1 H+ 12 C=4.438 


11.0 


317 


4 He+ 12 C=4.438 


3.00 


393 


1 H+ 14 N=4.438 


21.7 


291 


4 He+ 14 N=4.438 


9.64 


333 


1 H+ 16 0=4.438 


22.5 


156 


4 He+ ie O=4.438 


25.0 


223 


1 H+ 16 0=6.129 


13.0 


163 


4 He+ 16 0=6.129 


3.50 


269 


1 H+ 20 Ne=6.129 


19.0 


108 



(at energy -Epeak and with cross section <T pea k), in terms 
of the values of the reference cross section at the peak, 
o'peak; and at ni S h energy, o£,, given above: 



<j{E) = tTpeak 



^peak/ V-Epoak/ O, 



peak 



(13) 



For completeness, we give the values of -E pC ak and tT pe ak 
for the various cross sections considered in this paper in 
Table [TJ. 

In general, the error on the excitation cross sections 
and thus on the gamma-ray yields is typically 10% when- 
ever actual experimental data exist (this is the case for 
the two main lines of 12 C and 16 0), and of the order of 
20% to 50% when the values are simply estimated or ex- 
trapolated. This is quite substantial, especially when our 
goal is to look at line ratios, in the hope to determine 
the composition of the EPs and/or the ambient medium 
from gamma-ray line measurements. These errors, unfor- 
tunately, cannot be lowered but by increasing the exper- 
imental effort at terrestrial accelerators. This is strongly 
recommended in order to make the most of the opening 
field of gamma-ray astronomy. 

5.2. Photon yields for the 12 C, 14 N, 16 O, 20 Ne and 56 Fe 
"/-ray lines. 

In Fig. |a, we show the absolute 7-ray yields, A/< 3 - ;7 , cor- 
responding to the main 16 line at 6.13 MeV, for various 
projectiles and targets. 



The evolution of A/i J;7 as the injection energy increases 
can be interpreted in the following way. Photon produc- 
tion begins when E ln becomes greater than the reaction 
threshold. Then it increases sharply as E- m passes through 
the peak of the cross section, and increases more smoothly 
afterwards. As long as particle destruction or escape can 
be neglected, Eq. ( |l2| ) makes it clear that the number of 
photons produced is an increasing function of E- m , the 
upper bound of the integral. Physically, the particle pro- 
duces 7-rays all the way as its energy goes down from E- m 
to below the reaction threshold. If it is injected at a higher 
energy, it will produce 7-rays for a longer time, integrating 
the cross section over a larger range of energy. 

But when E- ln increases further, there comes a time 
when the projectile has a large probability of being de- 
stroyed (through a nuclear reaction) or escaping from the 
region under study (in a thin target model), before its 
energy drops below the reaction threshold. In this case, 
the effective energy range over which the cross section is 
integrated is reduced from below, and the overall 7-ray 
yield starts to decrease. For large enough E- ln , the particle 
never reaches the most efficient energy range correspond- 
ing to the peak of the cross section. Since both the destruc- 
tion and the excitation cross sections are roughly constant 
at high energy, the photon yield tends to an 'asymptotic 
value' where increasing the injection energy only shifts 
upwards the energy range of activity of the EP but does 
not change the integrated photon yield. This asymptotic 
value merely depends on the ratio of the excitation and 
destruction cross sections. 

Typically, for the main 7-ray lines to be expected in 
the ISM, one can see from Figs. || to ^| that several percent 
to up to 30% of the projectiles injected will produce a 
gamma-ray, with a peak of this number in the range E ln = 
100-300 MeV/n. 

While the decrease of Afi,j ;l at high energy is not very 
steep, it should be realized that the 7-ray production ef- 
ficiency, defined as the number of photons produced per 
erg of projectiles injected, is falling down more quickly, as 
shown on the right sides of Figs. || to |[ for the same reac- 
tions as on the left sides. The corresponding curves give a 
visual representation of the most efficient energy range for 
an EP to produce a given 7-ray line. They can be thought 
of as simple phenomenological tools: a simple look at them 
gives an idea of the kind of source spectrum and compo- 
sition required to reproduce any 7-ray line observational 
data. Note that contrary to what might have been naively 
expected, this range starts at an energy higher than the 
cross section peak, and extends to even higher energies. 
In other words, the most efficient way to produce gamma- 
rays in the ISM is to use EPs with energies between, say, 
10 to 300 MeV/n. 
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Fig. 9. Gamma-ray yields of a 1 H (left) and a 4 He (right) nucleus injected in a medium of solar metallicity, as a 
function of the injection energy. Contributions to various 7-ray lines are shown. 



5.3. Gamma-ray line emission synthesis 

From a practical point of view, the quantities A/ij; 7 (-Ei n ) 
also allow one to straightforwardly calculate the 7-ray line 
emission in a given astrophysical situation. It suffices to 
sum the contributions of each reaction involved, weighted 
according to the desired chemical abundances of both the 
source and the target. In other words, one can calculate 
the 7-ray line emission rate for any EP spectrum and 
composition in any medium (except maybe the most ex- 
tremely metal-rich), without needing to worry about par- 
ticle propagation and energy losses at all, as intended. 

Such a weighting is illustrated in Figs. || and [Io| where 
we show the 7-ray yields of H, He, C and O nuclei in a 
medium of solar metallicity. Note that although the 7-ray 
yields of C and O projectiles appear much higher than 
those of He (or H), they still have to be weighted by the 
relative abundances of the various projectiles among the 
EPs. The number of gamma-rays in the C and O lines 
produced by one EP injected in the ISM at energies above 
a few tens of MeV/n is typically between 10 -5 and 10 -4 . 

An interesting result is the fact that, in addition to 
12 C nuclei, 16 nuclei are also rather efficient in producing 
the 12 C line at 4.44 MeV. In Fig. [n], we have also shown 
the 12 C*/ 16 0* emission line ratio for the three projectiles 
producing both of these lines, namely H, He and O. This 
can be used to estimate quickly the probable composition 
of EPs producing any observed 12 C*/ 16 0* line ratio. 

6. Analytical estimates 

In an earlier work, Bykov and Bloemen (1994) calculated 
the average photon yield of a single energetic nucleus suf- 
fering from Coulombian energy losses (neglecting nuclear 
destruction). It is worth comparing our results with their 
analytical approximation, obtained with the Bethe-Bloch 
formula for energy losses, and extending their approached 



1 cr 



10 3 




10" 4 

10° 

Injection energy (MeV/n) 
Fig. 10. Same as Fig. O for 12 C and 16 projectiles. 



formula to situations where the excitation cross section 
follows a simple analytical law (such as that mentioned in 
Sect. [5.l| ), and where the nuclear destruction is important. 

Depending on the EP injection energy, one can identify 
two opposite 'regimes' where the estimation of the pho- 
ton yields A/ij; 7 (-Ein) can be simplified. From the physical 
point of view, this depends whether the particle energy 
losses can be neglected with respect to nuclear destruc- 
tion, or vice versa. On Fig. [l2| we have drawn the de- 
struction cross section, o~d {E) , together with the energy 
loss cross section, o~i oss (E), defined by: 



floss — 



1 



n v(E)Ti oss (E) 



\E(E)\ 
nvE 



(14) 
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Fig. 11. 12 c*/ 16 0* emission line ratio produced by H, He 
and O nuclei injected into the ISM with solar metallicity, 
as a function of the injection energy. 
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Fig. 12. Comparison of the total inelastic cross section 
(nuclear destruction) and the energy loss cross section for 
12 C and 1 H. The 12 C excitation cross section is also shown 
for comparison. 

for energetic 1 H and 12 C nuclei. Note that with this defi- 
nition, Eq. ( [l2|) can be rewritten in a simple way as: 



O"lo 



(15) 



It can be seen on Fig. [T^ that at low energy, destruction 
is negligible compared to energy losses, while the opposite 
is true at high energy. The exact transition energy depends 
on the nucleus considered, and ranges between 200 and 
300 MeV/n. This is consistent with the results of Fig. 
showing the transition of the survival probability from 1 
to around this energy. 



6.1. Low energy limit 

In the case when particle destruction can be neglected, one 
can set the survival probability, Vi(E- m , E), equal to 1 in 
Eq. ( |l5|) (this is the case studied by Bykov and Bloemen, 
1994). The Bethe-Bloch formula for the energy losses of a 
nucleus (Z,A) gives: 



Z 2 3 , ... m e c 2 
moss{E) = -j X -cr T ln(A) X — — 



(16) 



where ctt is the Thomson cross section and In A is the 
usual Coulomb logarithm (see e.g. Lang, 1999). It is a 
slightly increasing function of energy (~ 7.3 at 10 MeV/n, 
and - 12 at 1 GeV/n). 

Reporting in Eq. ( [To] ) and assuming that the excita- 
tion cross section is constant, one obtains the result of 
Bykov and Bloemen (1994) in the non-relativistic limit 
(E ~ m p v 2 /2): 



A/" 7 (£ in ) = 



A m p <7p ro d 
Z 2 6m c In A ctt 



(17) 



where the cross section can then be considered as a func- 
tion of energy, although this amounts to moving a pr od(^) 
outside the integral in Eq. (|l5|), which is of course im- 
proper in principle. 

The above formula can be improved by using an 
approached expression for the nuclear excitation cross- 
section, which allows an analytical integration of Eq. [l^ 
(with Vi(Ei n ,E) = 1). By approximating the cross sec- 
tions as <Jp m d(E) = a x E~ x + b above the peak (the 
values for a and b being deduced from Eq. (|13| ) and Ta- 
ble 0) , with a linear connection from the threshold energy, 
E t h, to the peak energy, £peak, one obtains: 

for Eth < E < £peak, 



A^in) 



A 2(7p Ca k (Ei, 



E 



Ui 



2E h 



E 



i.h 



Z 2 9cr T In A m e c 2 m p c 2 E pcllk - E t h ' 
and for E > E peak , 



(18) 



A/" 7 (-Ein) = A/ 7 (£'poak) 



A/Z 2 
3ctt hi A 



4E 2 



m c c 2 m v c 2 V 2 



( aE~ x b 
■ + 2 



(19) 



where we have used the non-relativistic relation between 
v and E, which is justified for the energy range under 
consideration. 



6.2. High energy limit 

In the other limit, at high energy, the destruction time is 
so small compared to the energy loss time that the parti- 
cles will be destroyed before they have lost any significant 
amount of energy. We can thus assume that the gamma- 
ray emission occurs at a constant energy, namely the in- 
jection energy. The problem is solved straightforwardly 
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in terms of the production and destruction timescales (or 
cross sections): 



TP 
Tprod 



"prod 
CD 



(20) 



Indeed, if one particle is injected at time t = 0, 
the number of remaining particles at time t is N(t) = 
exp(— t/rrj), and the gamma-ray production rate is 
diV 7 /di = N(t)/T pmd = exp(-£/T D )/T prod . Integrating 
over time, from t = to oo, gives the above result. 

6.3. Accuracy of the analytical formula 

The approached formulae Eqs. (|l8|) and (|l^) give an ap- 
proximation of the photon yields at low energy, which 
may be noted A/^le^), while Eq. ( pp[ ) is accurate for 
the high energy limit, A/" 7i he(£')- One can then propose 
an approached formula valid in the whole energy range: 



N^E) = [A^LE^r 1 +A/" 7 , HE (£)~ 



(21) 



In Fig. [L3|, the various above approached formulae are 
compared with the results of the previous section for the 
reactions 12 C+ X H^- 12 C* and 16 0+ 4 He O*. As can 
be seen, the various approximations give reasonably good 
results in their respective energy range. However, it should 
be stressed that except for the high energy, the accuracy 
of our formulas depends essentially on the accuracy of the 
cross section modeling. Since the real cross sections do not 



follow the simple analytical expression used in Sect. 6.1 



we cannot expect the analytical photon yields calculated 
in this section to be accurate to more than a factor of two 
or so (as we could observe from the whole data set). This 
can be very problematic when considering gamma-line ra- 
tios. Therefore, we strongly recommend to use the results 
of Sect, ^instead, which are also uncertain, but only inso- 
far as the cross sections are not know. They thus provide 
the best estimates of individual gamma-ray yields given 
the present knowledge on the excitation cross sections. 

7. Summary 

In this paper, we have presented an easier way to cal- 
culate gamma-ray line emission from energetic particle 
interactions in the ISM. It is based on a simple mathe- 
matical transformation whose physical interpretation has 
been given and which allows one to work with the source 
spectrum of the EPs rather than the propagated spec- 
trum. Therefore, one does not need to worry about energy- 
dependent and nucleus-dependent energy losses of the par- 
ticles, nor about their nuclear destruction in-flight, as they 
are taken into account once and for all through the calcu- 
lation of absolute photon yields. The latter are the number 
of photons produced in each of the nuclear de-excitation 
lines by a given nucleus injected in the ISM at a given 
initial energy. These photon yields have been given here 



for various projectiles contributing to the main 12 C, 14 N, 
16 0, 20 Ne and 56 Fe de-excitation lines. Numerical tables 
and electronic versions of the results are available from 
the authors upon request. These photon yields are to be 
used instead of the nuclear excitation cross-sections, and 
might be thought of as 'effective cross-sections' taking into 
account the specific effects of particle propagation in the 
ISM. They can be used to calculate the 7-ray line emission 
induced by EPs with any spectrum and any composition 
in any medium with a metallicity lower than a few tens of 
the solar metallicity. 

In addition to simplifying the calculation of gamma- 
ray line emission, the individual EP gamma-ray yields also 
provide a direct, intuitive tool to analyze gamma-ray line 
data from a phenomenological point of view, and construct 
an EP source spectrum and composition which could re- 
produce the intensity of the gamma-ray emission and the 
various line ratios. The results presented here correspond 
to a thick target model, which is relevant to most astro- 
physical situations for EPs of energy lower than a few 
hundreds of MeV/n. However, the same formalism can be 
used to calculate the total EP photon yields in a target 
with any escape length, be it energy-dependent or not. 
Once these yields have been calculated once, they can be 
used in any situation with the same escape length, for any 
particle spectrum and any EP and target compositions. 
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